Non-linear Resistivity of a Two-Dimensional Electron Gas in a Magnetic Field 
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We develop a theory of nonlinear response to an electric field of a two-dimensional electron gas (2DEG) 
placed in a classically strong magnetic field. The latter leads to a non-linear current-voltage characteristic at a 
relatively weak electric field. The origin of the non-linearity is two-fold: the formation of a non-equilibrium 
electron distribution function, and the geometrical resonance in the inter-Landau-levels transitions rates. We 
find the dependence of the current-voltage characteristics on the electron relaxation rates in the 2DEG. 
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A magnetic field applied to a two-dimensional electron gas 
(2DEG) changes the energy spectrum and dynamics of elec- 
trons. It leads to a modification of the transport character- 
istics of the 2DEG even at relatively weak magnetic fields, at 
which the Landau levels 1 1] are not resolved yet and the Quan- 
tum Hall effect |2] is not developed. The most well-studied 
modification of that kind is the Shubnikov - de Haas (SdH) 
oscillations in the resistivity of 2DEG. Its observation, how- 
ever, is restricted to fairly low temperatures, T < hcj^, so that 
the thermal broadening T of the electron distribution is small 
compared to the Landau quantization energy Swc. In the low- 
temperature limit, the strength of the oscillations is controlled 
by the Dingle factor, A - [-;7r/(wc'ro)]; it yields information 
about the "quantum" lifetime tq of the electron |fT] due to scat- 
tering off disorder. 

Recently it was realized ifstl that the effect of a magnetic 
field on the dc non-linear transport, unlike the SdH oscil- 
lations of the linear resistivity, is not confined to low tem- 
peratures. Oscillations of the differential resistivity with the 
magnetic field at a finite level of current observed in Ref. fstl 
persisted to quite high temperatures (about 4K), while the 
conventional SdH effect was fully smeared out by tempera- 
ture. This finding was confirmed in later experiments 10, 110] 
where the differential resistivity was measured both as a func- 
tion of the applied magnetic field and transport current. The 
oscillations of the nonlinear resistance were associated with 
the geometrical resonance in the electron transitions between 
the Landau levels ^ that arises from the commensurability 
of the period in the spatial oscillations of the density of states 
(DOS) and the diameter IR^ of an electron cyclotron trajec- 
tory. Although this was a plausible explanation of the effect, 
it remained unclear, why the oscillations are so weakly sen- 
sitive to the temperature and which parameters of the 2DEG 
control the amplitude of the oscillations. 

Another notable effect of magnetic fields on the nonlinear 
transport in 2DEG was reported in Ref. |5] and deals with the 
region of relatively small cuiTent densities. In that regime, a 
sharp drop in the differential resistivity was observed. The ef- 
fect was attributed to the modification of the electron energy 
distribution caused by the current 131 • Clearly, this modifi- 
cation depends on the energy relaxation rate, and the mech- 



anisms behind the observations reported in Refs. ^ and ^ 
seem quite different from each other. 

The goal of our work is to show that the two seemingly 
different phenomena are essentially two manifestations of the 
electron kinetics described by a standard Boltzmann equation 
for a weakly disordered 2DEG in the presence of electric and 
magnetic fields. We demonstrate that the low-current nonlin- 
earity {3\ is the consequence of the variation of the occupation 
factors of electron states: the non-equilibrium population of 
states renders the transitions normally contributing to the dis- 
sipative current ineffective. At high currents, the effect of the 
electric and magnetic fields on the electron motion becomes 
important. The oscillations in the I-V characteristic are asso- 
ciated with the geometric resonance in the electron transitions. 

We evaluate the dissipative component of the electric cur- 
rent density in a 2DEG placed in a perpendicular magnetic 
field Z? as a function of electric field characterized by the di- 
mensionless parameter 
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proportional to the ratio of the work of electric field associated 
with the displacement of the guiding center of a cyclotron tra- 
jectory by 27?c to the Landau quantization energy ticoc- The 
displacement occurs due to the electron scattering off an im- 
purity, and does not exceed IR^ in a single scattering act. This 
geometrical constraint leads to the oscillations of the cuiTent 
with ^; each oscillation corresponds to an increase by fico^ of 
the maximal energy acquired by an electron from the elec- 
tric field in a single scattering event. The maximal displace- 
ment of the guiding center is reached for the scattering angle 
71 (backscattering), thus the amplitude of oscillations is pro- 
portional to the corresponding scattering rate 1/T(:7r). 

The "prefeiTed" values of the energy absorbed by an elec- 
tron from the electric field are multiples of tico^ because of 
the oscillations of the electron DOS associated with the Lan- 
dau quantization. It is interesting to note however, that even a 
strong electric field does not result in developing a substan- 
tial modulation in the electron energy distribution with the 
period htj^- The reason for that is the dual role electric field 
plays. On one hand, it promotes the build-up of electron dis- 



2 



tribution at the energies corresponding to the maxima of the 
DOS. On the other hand, it increases the electron diffusion 
in energy space, the corresponding coefficient of the spectral 
diffusion being proportional to the Joule losses. The latter ef- 
fect wins over the former one, and the electron distribution in 
energy gets smoother at higher fields. As the result, the os- 
cillatory part of the I-V characteristic reflects the modulation 
of the electron transition rates with the field, rather the mod- 
ifications in the electron distribution function; the amplitude 
of oscillations provides information about the backscattering 
rate l/T{n), hardly accessible in other experiments. 

In this work we express the dissipative current in terms of 
the inelastic relaxation rate 1 / Ti„ and harmonics 1 /t„ of the 
elastic electron scattering rate 1 /t(6') on angle 6: 
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Typically, the "quantum scattering time" tq is short, Tq <si 
Tin. However, the transport relaxation time, defined as I/th = 
^/to-^/ti, may be in an arbitrary relation with Ti,,. We show 
that the measurements of the dissipative current as a function 
of ^ at small < yfro/Tia) and large (> 1) values of ^ reveal 
the rates of inelastic relaxation and of the back-scattering off 
disorder, respectively. 

In the following, we consider the limit of high temperatures 
T > Tidj^llip-, when the non-linear resistivity is observed, but 
the SdH oscillations are already suppressed We also limit 
our analysis to the case of "classically strong" magnetic fields, 
i.e., we assume that (jj^tq < 1 while (jL>cTa ^ 1- The former 
condition allows us to keep only the first harmonic in the os- 
ciUations of the DOS 

v(e) = vo ( 1 _ 2icos ^ I , e'"''^'^". (3) 

\ rtWc / 

The condition t^Itq » 1 is routinely met in semiconductor 
heterostructures, and the domain of magnetic fields 1 /rtr <s; 
Wc ^ 1 /■!"() is quite wide. 

The dissipative part of the electric current 
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is determined by the stationary electron distribution function 
f{s, (fi), which is the solution of the following kinetic equation 

- cJcd^fie, if) = St^{/(e, ^)) + Sti„{/(e, if)}. (5) 

Here (p is the angle the electron momentum makes with the 
direction of the electric field. The first term in the right hand 
side of Eq. (|5]i is the collision integral for electron scattering 
off disorder: 
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Here W^p^' - eERi^[smifi' - sin^] is the work of the electric 
field in the course of the shift R^z x [n^- - n^] of the guiding 




FIG. 1: Electron scattering off impurity changes the momentum di- 
rection from to n^i and the position of the guiding center shifts 
by R^z X [n^, - n^]. 



center of the cyclotron trajectory, see Fig. [TJ unit vector z is 
perpendicular to the 2DEG plane, and - {cos (fi, sin (p] 
is directed along the electron momentum. The rate of such 
scattering events is given by 1 / T((fi' - (fi) and characterized 
by its harmonics 1/t„, see Eq. (|2]). The imbalance between 
scattering "in" and "out" terms is determined by the difference 
of the corresponding distribution functions /(e + W^^i , (fi') and 

The inelastic relaxation at sufficiently low temperatures is 
dominated by the electron-electron interaction, represented by 
the inelastic collision integral 

Sti„{/(e)) = J de' J dEM(E,s,e') (7) 
X [f(e)f(s^)f(s')f(s'_) - f{s)f(e^)f(s')f(s'_)] . 

Here /(e) = 1 - /(e), e+ - s + E, e'_ - s' - E and 
M(E, s, s') describes the dependence of the matrix element 
of the screened Coulomb interaction on the transferred energy 
E and the electron energies e and e' . 

To the first order in 1 /Wc^a <k 1 we look for a solution of 
the kinetic equation Eq. (|5]) in the form: 
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The first term in Eq. ([S) is the result for A -Q, when the effect 
of quantization by magnetic field is neglected. The second 
term in Eq. (jSJ is proportional to A and oscillates with e/Wc- 
Here we assume that T » oj^, so that we can separate fast 
oscillatory dependence on energy of /(e, (fi) with period Tud^ 
and smooth energy dependence of fjie) on the scale of tem- 
perature T of the 2DEG. We assume that tiu^^ <k T and that 
/7-(e) is close to the Fermi distribution function. 

The oscillations amplitude /] (e) of isotropic in momentum 
space component of the distribution function Eq. ([8]) is deter- 
mined by 
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For the amplitudes A 12 of the anisotropic component in 
Eq. ^ we have 

^A, cos ^ = (sin^ + KM) ' (9b) 



^^A..-(sin,cos^^^^^^7C^,(.))/^^ (9c) 

Here (...) stands for averaging over angle variables tp and (p'. 
The kernels 'K^^^}\s) are given by 



^M) , X sin[27r(g + W^^O/^f^c] - sin[27re/^Wc] 
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In the left hand side of Eq. ( |9al i we used the linearized form of 
the inelastic collision integral Eq. O. The inelastic relaxation 
rate l/rin(e) of the oscillating component of the distribution 
function Eq. ^ was calculated in |j7|]: 
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where E-p is the Fermi energy and k - Aire^vo. Solving 
Eq. ( l9ab with respect to /] , we obtain 
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Next we substitute /i(e) into Eqs. ( |9b] i and ( |9cl ) and find 

Ai(£) = -(2eERJoj,T^) [Yi(0 rUe)) + TziO] > 

A2(£) = (2e£/?JWeTtt)r,(^,Tin(e)). 

Here we introduced the following notations: 

ri(^,ri„) _ [dy(0/d(f Y2(0 _ d^yiO 
Tr_i+To'-r(0' ^tr df ' 



(12b) 
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and J„{0 are the Bessel functions. 

The isotropic in momentum and oscillatory in energy com- 
ponent of the distribution function f{s,(p) results in the 
ri(f, Tin) contribution to the amplitudes A 1,2 of anisotropic 
part of f{s, (fi) I?!. The second contribution, containing TziO^ 
is coming from the second term in the r.h.s. of Eq. ( l9b] i. This 
contribution was studied in Refs. and originates di- 

rectly from the effects of electric fields on the collision integral 
for scattering off disorder. 

Substituting the distribution function Eq. ([8j with oscillat- 
ing components Eq. ( fT2l i into Eq. Q and integrating over en- 
ergy, we obtain the dissipative current: 

jd = (ruE[l + 2A^F(0l F(0 = 2ri(^, fi„) + Y2(0, (14) 

where crp - e^R^vollTa is the Drude conductivity at large 
Hall angle, Wc^tr » 1. To perform the integration and 



simplify further analysis we replaced the inelastic relaxation 
rate 1/Tin(e) by a parameter 1/fin, which can be chosen as 
l/Tin ~ 1/Ti„(e= r). 

Equation (fT4l ) for the dissipative current together with 
Eqs. ( I12al ) and (fT3l l for functions y(f), T\{^,Tin) and r2(^) 
constitute the central result of the paper. Below we discuss the 
properties of functions ri(^, fi„) and r2(^), which determine 
the non-linear response of the dissipative current j^, Eq. ( fT4b . 
Then we consider a specific model for 1 /t„ in Eq. (|2]i and an- 
alyze the non-linear behavior of the current within the model. 

For weak electric fields, ^ <K 1, we obtain the following 
expressions for functions ri(^, fin) and r2(^) expanding the 
Bessel functions to lowest order in ^^: 
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The term in the denominator of Eq. JlSal i is legitimate in 
the limit Tin ^ t^, which may take place at sufficiently low 
electron temperatures. We also note that Eq. ( I15al ) coincides 
with the result of Ref. [TJ in the absence of microwave fields. 
In the strong-field hmit, ^ » 1 , we find 

ri(f,fin) To cos^^ TiiO 1 sin 2^ 

oc — . oc — . (10) 
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where 1/T(7r) — Yju^'"" /^n is the back-scattering rate off 
disorder. The r2-contribution, arising from the effect of 
electric field on the collision integral, is larger than the 
Fi -contribution, which arises from the stationary out-of- 
equilibrium component of the distribution function. The lat- 
ter contribution not only decays faster with the increase of ^ 
than the former one, but also contains an additional parameter 
to/t(7t) which is small for smooth disorder The amplitude of 
oscillations of current Eq. (fT4l i decays proportionally to 1 at 
^ » 1, but the oscillations in differential conductivity 

cr = djd/dE oc d[^r2(0]/di oc COS 2( 

do not vanish; its maxima and minima are situated at ^ = 
7Tk/2 with integer k. At smaller values of term TiiO ^Iso 
contributes to cr and results in dependence of the oscillations 
amplitude on ^ as well as in some shift of maxima and minima 
from ^ = 7rk/2. 

To discuss the properties of the non-linear current Eq. ( fT4l ) 
in a broad range of electric fields, we consider a specific model 
for the harmonics 1/t„ of the scattering rate off disordered 
potential due to charged impurities inside or in the proximity 
of2DEG: 



1 1 1 6„fi 
— = -; 2+ ' /^^l- 

T„ Tsm 1 +Xn'- Tsh 



(17) 



Here I/tsh is the scattering rate off impurities inside the 
2DEG, which produce sharp (5-correlated) potential for elec- 
trons. Charged impurities in the proximity of 2DEG produce 
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FIG. 2: Function F(() for disorder described by Eq. l ll7b with 
Tsm = Tsh/SO for several values of fi„ and^: fi„ = IOtsi,,^ = 0.02 
(solid); Tin = 10Tsh,;r = 0.01 (dashed); f,„ = 0.5T,h,x = 0.01 (dot- 
ted). Note that F(^) has a peak at small ( determined by the inelastic 
relaxation time f and oscillates at large Zf. At intermediate ^ ~ 1 
the positions of maxima and minima of F(f ) depend on the relation 
between various parameters of the model. 



oscillatory behavior of F(^). 

We calculated the dissipative component of electric current 
in response to the applied dc electric field of arbitrary strength 
within self-consistent Born approximation [1]. We show that 
the non-linear component of the current consists of two con- 
tributions. One contribution arises due to the formation of 
the out-of-equilibrium component of the distribution function, 
oscillating as a function of energy. The second contribution 
is the result of modification by electric field of electron scat- 
tering amplitudes off the disorder potential. We showed that 
the former contribution is important at relatively weak fields, 
while the latter one dominates in the high-field domain. There, 
the non-linear contribution to the current oscillates as a func- 
tion of the applied electric field. The amplitude of oscillations 
of the differential conductivity does not decrease with the in- 
crease of electric field (and at fixed magnetic field). It may 
be necessary to take into account the effect of heating on the 
quantum scattering time and thus on the Dingle factor A in or- 
der to explain the suppression of oscillations observed in j^. 



a smooth potential resulting in electron scattering on small an- 
gle 9 ~ yjx <K 1, where ;t' can be estimated x ^ 
with Ap and ^ being the Fermi wavelength and the correlation 
length of the disorder potential, respectively. The two restric- 
tions for validity of Eq. ( fT4l i on the strength of magnetic field 
(WcTti » 1 and cl>cTo < 1) can be satisfied simultaneously for 

For the disorder characterized by harmonics of the scat- 
tering rate Eq. (flTl l. we obtain y(f) = JoiO/'''sh + 
l/C^sm Vl +X(^)- Here we omitted the term, which arises 
from the back-scattering off" smooth disorder and yields ex- 
ponentially small (exp(-7T/ sjx)) contribution to 1/T(7r), cf. 
Eq. ( fTSI l. Substituting this expression for 'x(^) in Eq. iT3[ . we 
can evaluate the current Eq. (fT4l) at arbitrary At ^ <K 1 / sjx 
we have 

Tu- fsn, fr'+T;,'[l-y2(^)] + f-J,f2/2 ' 

(18) 

where l/fsm - A'/^sm is the smooth disorder contribution 
to the transport scattering rate. Equation ( fTSl l covers both 
the regime of relatively weak fields, where inelastic scatter- 
ing is important, and the regime of strong fields, exhibiting 
prominent oscillations. Function F(^) has a sharp feature at 
^ ~ V'!"sh/'!'in- At these fields, the spectral diffusion of elec- 
trons caused by electric field becomes comparable with the 
inelastic relaxation. At stronger fields, ( > V^smAHn^ the two 
kinds of disorder, smooth and "sharp", yield two separate con- 
tributions to F(0, 

F(o _ 1 - - 2(1 + yrr"^) [JKor 

Here only the sharp component of disorder contributes to the 



Finally, we considered the limit oJcTq < 1, and therefore 
assumed that the oscillations in the DOS are described by 
one harmonic with period in energy Swc, Eq. (O. In stronger 
fields, the DOS remains periodic in energy^ with the same pe- 
riod hciJc, but contains higher harmonics llK^. These higher 
harmonics in the DOS result in a more complicated form of 
the oscillatory part of the non-linear resistivity. 
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